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A brief report on the QCD factorization approach to exclusive B decays and its applications is presented. In 
an appendix a subtle issue concerning partial integration in the parameterization of infrared-sensitive power 
corrections is clarified. 



1. INTRODUCTION 

Exclusive non-leptonic decays play a crucial 
part in the on-going programme to clarify the 
sources of CP violation and rare flavour-changing 
phenomena. As regards CP violation this is illus- 
trated by the fact that there exists a fair knowl- 
edge of the length of each of the sides of the uni- 
tarity triangle from semi-leptonic B decays and 
BB mixing. Since recently, the angle [3 (more 
precisely, sin 2(3) is determined fairly accurately 
by the time-dependent CP asymmetry in the non- 
leptonic decay Bd — ► J/ipKs and related ones. 
However, the direct measurement of the other an- 
gle 7, the phase of V* b in the standard phase con- 
vention, remains a challenge to theory and ex- 
periment. Although 7 is determined indirectly 
through the other measurements, the complemen- 
tary direct determination is important, since me- 
son mixing and rare non-leptonic decays could be 
affected by new flavour-changing interactions. 

The angle 7 can be determined from decays 
with some kind of interference of b — > cuD (no 
phase) and b — > ucD (phase 7) tree transitions 
and their conjugates (D — d, s) or from decays 
with interference of b — > uuD (tree, phase 7) and 
b — > Dqq (penguin, phase (D = s), /3 (D = d)). 
There exist a few methods that allow one to ob- 
tain 7 in a theoretically clean way, which means 
that no calculation of strong interaction effects is 
needed, but these methods are all difficult to real- 
ize experimentally at present. The more general 
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situation is that some information on the decay 
amplitudes is required to extract 7 from the inter- 
ference term in branching fractions or CP asym- 
metries. In technical terms one needs to compute 
the matrix elements (f\Oi\B), where Oi is an op- 
erator in the weak effective Hamiltonian and / is 
a particular exclusive final state. 

2. THEORY OF EXCLUSIVE DECAYS 

2.1. QCD-improved factorization 

The QCD factorization approach 0J^] uses 
heavy quark expansion methods (7715 3> Aqcd) 
and soft-collinear factorization (particle energies 
^ Aqcd , "colour-transparency" ) to compute the 
matrix elements (f\Oi\B) in an expansion in 
l/m& and a s . Only the leading term in I/7715 
assumes a simple form. The basic formula is 

{M 1 M 2 \O l \B) = F B ^ Mi (0) f duT'{u)^ M M 

Jo 

+ [dtdudvT II (Z > U,v)$ B (0$M 1 (v)$M 2 (u), (1) 

valid whenever M 2 is a light meson (or, at least in 
principle, a charmonium state) and M\ any light 
or heavy meson. (M\ is the meson that picks up 
the spectator quark from the B meson. If Mi is 
heavy the second line is higher order in 1/7715.) 
The formula shows that there is no long-distance 
interaction between the constituents of the me- 
son M2 and the (BM%) system at leading order 
in I/7715. This is the precise meaning of factor- 
ization. In particular, the hard scattering kernels 
T 1 - 11 contain all rescattering phases at this or- 
der and include hard interactions with the specta- 
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tor quark (T 11 ). The following non-perturbative 
inputs are needed: F B ^ Ml - a B — > M\ form 
factor; <j>M; - the standard light-cone distribu- 
tion amplitudes for light mesons; 3>b - the light- 
cone distribution amplitude of the B meson ||- 
, some aspects of which remain to be under- 
stood. The QCD factorization approach extends 
the Brodsky-Lepage approach to exclusive hard 
processes [0, because it shows factorization also 
when soft interactions dominate the B — > M\ 
form factor, as is expected ||. 

2.2. Theory status 

Leading order in 1/mb- Proofs of the factor- 
ization formula for heavy-light final states exist 
at two loops ||] and probably to all orders PJlO| . 
These proofs must be extended to charmless fi- 
nal states by including the hard-spectator inter- 
action. In practical applications the calculations 
are currently done at next-to-leading order in per- 
turbation theory. There exist theoretical uncer- 
tainties from input parameters, which can be esti- 
mated, the most important of which are \ V u b\ and 
form factors. The theory is presently complete for 
flavour non-singlet mesons. 

Power corrections in l/m&. Factorization does 
not and is not expected to hold at subleading 
order. Attempts to compute subleading power 
corrections to hard spectator-scattering in per- 
turbation theory usually result in infrared diver- 
gences, which signal the breakdown of factoriza- 
tion. A subset of power corrections has been 
parametrized or estimated in ||Jll|,|l^] . It is com- 
mon practice in applications of perturbative QCD 
to not include errors from power corrections. One 
would wish to do better here, first because non- 
leptonic decays test fundamental aspects of weak 
interactions, second because some power correc- 
tions related to scalar currents are enhanced by 
factors such as m^/((m„+m d )AQCD) @- At least 
these effects should be estimated and included 
into the error budget. This has been attempted in 
fj^ ]. All weak annihilation contributions belong 
to this class of effects. 

2.3. Overview of applications 

QCD-improved factorization has been applied 
to date to the following final states: 



- PP, where P denotes a light pseudo-scalar 
non-singlet meson |l|Jl2]-[r^] . This is the ap- 
plication developed furthest. 

- some PV final states, where V denotes a 
light vector meson |l^ , p^ ]. 

- some VV final states (In the heavy 
quark limit both vector mesons are longitu- 
dinally polarized.) 

- D(*'L, where L denotes a light pseudo- 
scalar or vector meson P. p^pO] . 

- D(*'X, where X denotes a scalar or "ex- 
otic" meson 

- J/ipK |§,|§, J/ipK* |4|||. The factor- 
ization approach should apply in the limit 
m c — > oo, but one expects corrections to 
factorization of order Aqcd/(™c«s) ~ 1 ]§]■ 

- Radiative final states V-/ [§|-§i and Vl + l~ 
p6|] including isospin-breaking effects that 
could cause a different rate for neutral and 
charged B meson decay |29f| . 

- New physics effects due to a modification of 
the b — > sg transition |30| or due to R-parity 
violating supersymmetry pl| ] . 

Of these final states only D^L and PP will be 
discussed below. 

2.4. Other approaches 

Naive factorization p2[ . Naive factorization 
follows from the factorization formula for D^L, 
L\L,2 final states, when a s — > and m& — * oo. 
QCD-improved factorization extends naive fac- 
torization in the same sense in which the QCD- 
improved parton model extends the naive parton 
model. 

The "PQCD" approach Q. In this approach 
the decay amplitude is calculated in the Brodsky- 
Lepage formalism. The main difference with 
the QCD-improved factorization approach arises 
through the stronger assumption that Sudakov 
suppression renders the B — > M\ form factor and 
many power-suppressed effects (such as weak an- 
nihilation) calculable. The advocates of the QCD 
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factorization approach do not agree with the pro- 
ponents of the PQCD approach on conceptual as- 
pects of the implementation of Sudakov suppres- 
sion and on the claim that the result is perturba- 
tive and accurate |Q. 

3. SELECTED RESULTS 

3.1. B -> £>(*>£, 

The phenomenology of these decays is dis- 
cussed in [||. QCD-improved factorization ex- 
plains the experimentally well-established fact 
that the colour-allowed class-I decays to D^ + L~ 
are controlled by an approximately universal coef- 
ficient | ai| ~ 1.05. On the other hand, factoriza- 
tion is not expected to hold for colour-suppressed 
class-II decays to D^*'°L°. The corresponding co- 
efficient ai is therefore depends on the final state. 
The precise nature of ai depends, however, on the 
counting scheme. If the D meson is considered as 
heavy (m c /mb fixed, nib — * oo), xa2ja\ is power- 
suppressed in the heavy quark limit. The power 
suppression is due to 

T _ {m%-ml)f D F ^(ml) 

but with x ~ 0.9 this suppression is only formal. 
If the D meson is considered as light (m c fixed, 
TOfc — > oo), xaila\ counts as order 1 and becomes 
calculable. (A crude estimate has been given in 
j2|.) In this case one expects an important non- 
universal contribution from spectator-scattering. 
In both cases xa^/ax should be smaller than 1 
due to colour factors. 
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Figure 1. Determination of xaz/ai for D^n 
final states |^0| . 



The measurement of B° — > Z}(*)°7r° allows one 
for the first time to construct the complete isospin 
triangle for B — > D^tt decays and hence to de- 
termine the magnitude and rescattering phase of 
xa-i /a±, see Figure [l| The result is that |«2 1 — 0.4 
is sizeable for ttD final states. Together with 
|o2 1 — 0.25 for KJ/ip final states this provides 
clear evidence for the non-universality of colour- 
suppressed decays, but it is difficult to interpret 
this result quantitatively since the theory is not 
expected to be accurate for ir°D° and KJ/ip final 
states. For further discussion of this result from 
the perspective of QCD factorization see [j2(| . 

3.2. 7 from B — * tttt,ttK 

The possibility to determine the CP-violating 
angle 7 by comparing the calculation of branch- 
ing fractions into 7777 and irK final states with the 
corresponding data has been investigated in de- 
tail j[2| (see also The branching fractions 
for the modes B T -> tt^tt and B T -> ir^K , 
which depend only on a single weak phase to very 
good approximation, are well described by the 
theory. This demonstrates that the magnitude 
of the tree and penguin amplitude is obtained 
correctly, where for the penguin amplitude the 
1-loop radiative correction is important to reach 
this conclusion. There is, however, a relatively 
large normalization uncertainty for the irK final 
states, which are sensitive to weak annihilation 
and the strange quark mass through the scalar 
penguin amplitude. This uncertainty can be par- 
tially eliminated by taking ratios of branching 
fractions. The agreement is less good for branch- 
ing fractions with significant interference of tree 
and penguin amplitudes, if 7 is assumed to take 
values around 55° as favoured by indirect con- 
straints. This may indicate an interesting alter- 
native determination of 7, but the CP-averaged 
branching fractions sensitive to 7 also depend on 
the strong phase difference of the two interfer- 
ing amplitudes. This dependence is suppressed 
as long as the phase difference is small since 
only its cosine enters the calculation. With some 
exceptions (such as the final state 7r°7r°) QCD- 
improved factorization predicts that strong phase 
differences are not large, since they are paramet- 
rically of order a s or AQco/^b, and it therefore 
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Figure 2. 95% (solid), 90% (dashed) and 68% 
(short-dashed) confidence level contours in the 
(p, fj) plane obtained from a global fit to the CP 
averaged B — > ttK, tttt branching fractions, us- 
ing the scanning method as described in p5| . 
The right dot shows the overall best fit, whereas 
the left dot indicates the best fit for the default 
hadronic parameter set. The light-shaded region 
indicates the region preferred by the standard 
global fit [Q, excluding the direct measurement 
of sin(2/3). 



also predicts small direct CP asymmetries, up to 
(10 — 15)% (with a preference around 5%) for the 
it® , ir^K final states. The current data seem 
to favour small CP asymmetries, but they are not 
yet accurate enough to decide upon whether the 
QCD factorization approach allows one to predict 
strong phases quantitatively. 

The current interpretation of charmless non- 
leptonic B decays can be summarized by a global 
fit of the Wolfenstein parameters p, fj to the six 
CP averaged 7T7t, nK branching fractions shown 
in Figure |[ The result is consistent with the 
standard fit based on meson mixing and \V u b\, 
but shows a preference for larger 7 or smaller 
\Vub\- If the estimate of the theory uncertainty 
(included in the curves in the Figure) is cor- 
rect, non-leptonic decays together with \ V u b\ from 
semi-leptonic decays already imply the existence 
of a CP-violating phase of V u b at the 2-3 a level. 

3.3. sin(2a) or 7 from Bd(t) — ► 7r+7r _ 

The angle a (or 7, given the BB mixing phase), 
can be determined from the time-dependent CP 
asymmetry in Bd — » 7r + 7r~ decay alone, if the rel- 
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Figure 3. Constraint from S V7r in the (p,fj) 
plane. 



ative magnitude of the penguin amplitude, P/T, 
can be computed. The asymmetry is given by 



A CP [irir}(t) = S^sin(AM Bd t) 

+ A^[mr] cos{AM Bd t), 



(3) 



where S n7r = sin(2a), if P/T — 0, in which case 
the direct CP asymmetry ^cp[7r7r] vanishes. The 
direct CP asymmetry is proportional to the sine 
of the strong phase of P/T and can be used as 
a phenomenological check of the computation of 
P/T . Figure || displays how a measurement of 
Sun constrains (p,fj), when P/T is computed in 
the QCD factorization approach |l2|. The Fig- 
ure illustrates that even if theoretical (or experi- 
mental) uncertainties prevent an accurate deter- 
mination of sin(2a) in this way, the inaccurate 
result on sin(2a) still translates into a useful con- 
straint in the (p, fj) plane. This reflects the fact 
that other observables do not currently constrain 
sin(2a) very well. 

4. CONCLUSION 

Interpreting the new data on exclusive two- 
body decays of B mesons in terms of fundamen- 
tal parameters of the Standard Model requires 
- at least for the present - control over strong 
interaction effects. QCD-improved factorization 
allows us to formulate and discuss the problem 
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systematically in the expansion parameters a s 
and Aqcd /mb- The first comparisons of the the- 
ory with data on heavy-light and charmless fi- 
nal states are encouraging, but the theory may 
not always be accurate or useful. More data and 
calculations should help understanding where the 
theory works or fails, and if it fails, why it does. 
So far no unexpected failures are observed. If 
the estimates of theoretical uncertainties are cor- 
rect, B — » tttt,ttK together with \V u b\ from semi- 
leptonic decay establishes 7 ^ at the 2-3 a level. 
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APPENDIX 

In this appendix I address a subtle point in 
the parameterization of infrared sensitive twist- 
3 corrections to non-leptonic decays. The point 
concerns the legitimacy of integrating by parts in 
convolution integrals with apparently divergent 
boundary terms. Depending on the procedure the 
result appears to be either a linear or a logarith- 
mic infrared divergence. The difference is crucial, 
since a linear divergence compensates the l/rrib- 
suppression and hence spoils factorization even 
in the heavy quark limit, while a logarithmic di- 
vergence does not. The results of the following 
discussion have been used but have not been ex- 
plained in our previous work |L2| . 

In order to dispose of the need to introduce 
B meson distribution amplitudes I consider an 
analogous problem that appears in the calcula- 
tion of the pion form factor {n + (p')\u , y tl u\TT + (p)} 
at large momentum transfer Q 2 — — (jp 1 — p) 2 . The 
light-cone expansion in 1/Q 2 can be implemented 
by multiplying the on-shcll four-quark scattering 
amplitude with a projection operator. For the fol- 
lowing discussion it will suffice to restrict the pro- 
jection to the two-particle quark-antiquark state 
and to neglect meson mass corrections. The co- 
ordinate space definitions [|6| 

{TT + (p')\u a (y)d (x)\Q) = I due ^vi.y+u P >. x) 



x { /is <j>(u) - ^75 ( <j>p(u) - cr piv p' fl z v ^" U 



0a 



are equivalent to projecting the scattering ampli- 
tude with 



0a 



^'75 <p(u) - (MrfB ( 4> P {u) 



, M (j> g(u) 8 

6 dk v 



0a 



and to integrate over u from to 1. (All quark- 
antiquark operators are assumed to be colour sin- 
glets.) I will refer to this projection as "PI" . Here 
k is the momentum of the quark in the pion and 
the antiquark momentum is p' — k. The deriva- 
tive can be decomposed into a transverse deriva- 
tive, a term proportional to p' v , which does not 
contribute, and a term proportional to p v . Us- 
ing u — p ■ k jp 1 ■ p one performs an integration by 
parts on this term and neglecting boundary terms 
(because 4> a {u) vanishes at and 1), one obtains 
the projection "P2" : 



M 0a 



if* 



■ IO av —. ; V ICT^p p 



p -p 
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6 dk±„ 
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The calculation of the form factor then gives 

2 r2 /•! 



(tt (p )|w7m«I 71 " (p)) 



9Q 2 



dudv 



4>{U)4>{V) , fJ^_ rp 



where the first term in the bracket represents the 
leading-twist result and T4 a twist-4 correction, 
which is the quantity of interest here. 

The explicit form of T4 depends on whether PI 
or P2 is used for the calculation. With PI the 
result is 



4> p {u) 



— 4> P (v) + — ^-t 1 } 



() 



+ 



=2 ^( V ) 



2_ <Mv) 

uv 3 6 



+ (u <-» v,p <-» p') 
while using P2 leads to 



4>p{u) 



— <t>p{v) + — + — >P» 

uv uv uv 



G 



+ 



+ 



1 1 , s 1 



2 ft f P M 



The two results are identical, when an integra- 
tion by parts is performed and when the bound- 
ary terms are zero. This does not seem to be 
justified, however, since the boundary terms are 
non-zero, in fact infinite, when the asymptotic 
distribution amplitudes <fi p (w) = 1, </v(w) = 6ww 
are inserted.^] Not performing the integration by 
parts, it is seen that the result from PI is lin- 
early divergent as v —> 1, while the result from 
P2 exhibits only a logarithmic divergence. Since 
non-perturbative infrared effects regulate the di- 
vergence at v A/Q (when the quark momen- 
tum becomes of order A), the Pl-result is of or- 
der Q/A, parametrically larger than the P2-result 
which is only of order In Q/A. This can be made 
more explicit by regulating the convolution in- 
tegrals through introducing a gluon mass A. The 
main effect of this is to replace uv by uv+£ where 
£ = X 2 /Q 2 . The convolution integrals can then 
be done with the two results: 

J dudvT 4 (Pl) = f-i-2Li 2 f-i^ (p+p% 

dudvT A (P2) = -2lA 2 (-~ \ (p + p')» 

The difference — 1/£ can be attributed to a lin- 
early divergent boundary term 

1 

T 



lim — — — — 

v^o v(uv + c,) 6 



Now it is clear that was has gone wrong here 
is that one divides by 0. The numerator is 
</> CT (l) = 0, which is always correct, but the de- 
nominator zero arises from the virtuality of a 
quark propagator, which approaches zero and 
which is not regularized by introducing a gluon 
mass. To make the calculation well-defined all 
small- virtuality regions must be properly regular- 
ized. One possibility is to assign a small mass to 

In the present approximation that neglects quark-anti- 
quark-gluon amplitudes the asymptotic distribution am- 
plitudes are exact and follow from the equation of motion 
constraints. A possible caveat here is that the solutions 
given above also involve an integration by parts in which 
boundary terms should be zero even when the equation of 
motion is inserted into any Green function. The following 
discussion shows that this assumption can be made. 



the quarks as well. Denoting the corresponding 
parameter by ry, this leads to the boundary term 

i <M«) 



lim 



0. 



i)->0 (v + Tj)(uV + £) 6 

It is not difficult to see that in this particular case 
the amplitude has a finite limit for rj — > (but £ 
still non-zero) after the convolution integrals are 
performed. The only effect of the complete regu- 
larization is therefore to make all integrations by 
parts well-defined and to eliminate the spurious 
linearly divergent boundary term. We can there- 
fore conclude that the correct result is 

T 4 (P2) = A(p + p ') 
uv 

which exhibits a double logarithmic divergence, 
one from each convolution integral. This coin- 
cides with the result for the chirally-enhanced 
twist-4 corrections obtained long ago by Geshken- 
bein and Terentev |57| . 

The discussion can be analogously repeated for 
non-leptonic B meson decays, where the problem 
appears in the computation of l/m^-suppressed 
terms. The conclusion is again that integration 
by parts is allowed and that the correct result 
is obtained directly with the projector P2. On 
the other hand using PI requires one to regular- 
ize carefully all propagators that can go close to 
mass-shell near the endpoints of the integration 
region. 

There is a deeper reason why the infrared di- 
vergence in the twist-4 correction to the pion form 
factor should not be power-like. There exist two 
basic scattering mechanisms that contribute to 
the form factor. One is the hard scattering mech- 
anism, in which all components of a given Fock 
state of the pion participate in a hard interaction. 
The other is a soft-overlap contribution, in which 
case a number of Fock components with small 
longitudinal momentum fraction does not partic- 
ipate in a hard interaction. A complete asymp- 
totic expansion of the form factor should com- 
bine both mechanisms consistently. Recent work 
on the analogous problem for the B — > 7r form 
factor [El together with the structure of the ra- 
diatively corrected QCD sum rule for the B — > it 
form factor |38| suggests the following factoriza- 
tion formula that incorporates the hard scattering 
and soft-overlap contribution: 

1 



■ T2 * (j>2 * ' 



7 



+ caq 2 , aO + r 4 rcg (/i) * 03 * 03 + o(i/Q 6 



Here 02 and 03 represent twist-2 and twist-3 
light-cone distribution amplitudes and £,(Q 2 ,/J,) 
a soft-overlap form factor which scales as 1/Q 4 
and which can be defined in an effective theory 
in which all hard scatterings are integrated out. 
Tl C6 (fi) denotes the hard-scattering kernel at the 
twist-4 level discussed above except that the end- 
point divergence is factorized into the soft-overlap 
form factor, such that both terms depend on a 
factorization parameter /i. If this (conjectured) 
factorization formula is to be correct, the end- 
point contribution of the twist-4 term can be at 
most logarithmically enhanced. 
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